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ABSTRACT 



Although the cubic "small universe" has been ruled out by COBE /DMR results 
as an interesting cosmological model, we still have the possibility of living in a universe 
with a more anisotropic topology such as a rectangular "small universe" with one or 
two of its dimensions significantly smaller than the present horizon (which we refer to 
as T^- and T^-models, respectively). In order to rule out these anisotropic topologies 
as well, we apply a new data analysis method that searches for the specific kind of 
symmetries that these models should produce. We find that the 4 year COBE/DMR 
data set a lower limit on the smallest cell size for T^- and T^-models of 3000/i~^Mpc, at 
95% confidence, for a scale invariant power spectrum (n=l). These results imply that 
all toroidal universes (cubes and rectangles) are ruled out as interesting cosmological 
models. 



Subject headings: cosmic microwave background, large-scale structure of universe. 



-3- 



1. INTRODUCTION 

In the past few years, mainly after the discovery of CMB anisotropies by COBE/DMR (Smoot 
et al. 1992), the study of the topology of the universe has become an important problem for 
cosmologists and some hypotheses, such as the "small universe" model (see e.g. Ellis and Schreiber 
1986), have received considerable attention. From the theoretical point of view, it is possible to have 
quantum creation of the universe with a multiply-connected topology (Zel'dovich and Starobinsky 
1984). Prom the observational side, this model has been used to explain the "observed" periodicity 
in the distributions of quasars (Fang and Sato 1985) and galaxies (Broadhurst et al. 1990). 

Almost all work on "small universes" has been limited to the case where the spatial sections 
form a rectangular basic cell with sides Lx,Ly, and with opposite faces topologically connected, a 
topology known as toroidal. The three-dimensional cubic torus is the simplest model among all 
possible multiply-connected topologies, in which all three sides have the same size L = = Ly = 
Lz- In spite of the fact that cubic T'^-model has been ruled out by COBE results (Sokolov 1993; 
Starobinsky 1993, hereafter S93; Stevens et al. 1993; Jing and Fang 1994; de Oliveira-Costa and 
Smoot 1995, hereafter dOCS95), the possibility that we live in a universe with a more anisotropic 
topology, such as a rectangular torus T^, is an open problem that has not been ruled out yet. 
For instance, if the toroidal model is not a cube, but a rectangle with sides / Ly / and 
with one or two of its dimensions significantly smaller than the horizon Rh (^ 2cHq^), this small 
rectangular universe cannot be completely excluded by any of the previous analyses: constraints 
from the DMR data merely require that at least one of the sides of the cell be larger than Rh. 

As pointed out by S93 and Fang (1993), if the rectangular T^-universe has one of the cell sizes 
smaller than the horizon and the other two cell sizes are of the order of or larger than the horizon 
(for instance, ^ Rh and Ly,Lz ^ Rh), the large scale CMB pattern shows the existence of a 
symmetry plane, and if it has two cell sizes smaller than the horizon and the third cell size is of the 
order of or larger than the horizon (for instance, Lx,Ly ^ Rh and > Rh), the CMB pattern 
shows the existence of a symmetry axis. We call the former case a T^-modcl because the spatial 
topology of the universe becomes just in the limit Ly, ^ oo with being fixed. The later 
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case is denoted a T^-model for the same reason (the corresponding Hmit is — >^ oo with L^,Ly 



being fixed). It is clear (and our calculations confirm it) that dependence of CMB fluctuations on 
any cell size is very small once it exceeds the horizon. In previous work (dOCS95), we computed 
the full covariance matrix for all multipole components and used a x^-technique to place a lower 
limit on the cell size L of the cubic T^-models. However, we cannot apply this same approach to 
study the T^- and T^-universes. As we explain in the next section, the observed power spectrum 
of these models depends not only on the cell size but also strongly on the cell orientation relative 
to the Galaxy cut. 

Our goal is to show that the COBE /DMR maps have the ability to discriminate and rule 
out T^- and T^-models. We use a different approach to study these models in which we constrain 
their sizes by looking for the symmetries that they would produce in the CMB, obtaining strong 
constraints from the 4 year COBE /DMR data. 



If the density fluctuations are adiabatic and the Universe is spatially flat, the Sachs- Wolfe 
fluctuations in the CMB are given by 



(Peebles, 1982), where r is the vector with length = 2cHq that points in the direction of 
observation {6,(p), Hq is the Hubble constant (written here as lOO/i km s~^ Mpc~^) and 6]^ is the 
density fluctuation in Fourier space, with the sum taken over all wave numbers k. Here we neglect 
the difference between the horizon surface and the surface of the last scattering, which is justified 
for / < 30. 

It is customary to expand the CMB fluctuations in spherical harmonics 



2. 



SYMMETRIES IN THE CMB DUE TO TOPOLOGY 




(1) 




1=0 m=-l 



(2) 
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where a;^ are the spherical harmonic coefficients and r is the unit vector in direction r. The 
coefficients aim are given by 



aim = -27Ti^^Y.^jl(kRH)Yi*m(k), (3) 
k 



where ji are spherical Bessel functions of order /. If we assume that the CMB fluctuations 6T/T 
are a Gaussian random field, the coefficients aim are Gaussian random variables with zero mean 
and covariance matrix 

'-^jl{kRH)jl'{kRH)Ylm{m'm'{^)- (4) 

In a Euclidean topology the universe is isotropic, the sum in (Q) is replaced by an integral and the 
power spectrum Ci is related to the coefficients aim by 

{alm^Vm') = Cl^U'^mm' (5) 

(see e.g. Bond and Efsthatiou 1987). However, in a toroidal universe this is not the case. In this 
model, only wave numbers that are harmonics of the cell size are allowed. We have a discrete k 
spectrum 

(6) 



_ ^ ( Px Py Pz 
Rh \Rx Ry Rz 



(Zel'dovich 1973; Fang and Houjun 1987), where px, Py and Pz are integers and Rx = Lx/Rh, 
Ry = Ly/Rn and R^ = L^/ Ru- 
in previous work (dOCS95), we set limits on the cubic T^-models assuming that, for a given 
cell size, the quantity Ci = giqrr E I^ZmP was fairly independent of the cell orientation, even with 
a 20° Galaxy cut. In other words, if Ci is almost independent of the cell orientation, we can make 
the approximation that all cell orientations for that given cell size can be simultaneously ruled out 
by a x^-t6st on the Ci coefficients and, in that way, test our model just considering changes in 
the cell size L. However, in the case of more strongly anisotropic cell configurations such as T^- 
and T^-models, the quantity Ci does depend on the cell orientation and the x^-test on the power 
spectrum cannot be used anymore. If we try to apply the power spectrum method to these models. 
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it will require testing a six parameter family of models with three parameters corresponding to the 
cell orientation in addition to the cell sizes L^, Ly and L^. 

In order to illustrate these anisotropic cell configurations, we plot a realization of two extreme 
cases: a T^-universe (Figure lA, upper left) with dimensions {Rx, Ry, Rz) = (3,3,0.3) and a T^- 
universe (Figure IB, upper right) with dimensions {Rx, Ry, Rz) = (0.3,0.3,3). Both models are 
plotted in galactic coordinates and have a scale invariant power spectrum (n=l). From equations 
dl]) and (^), we see that when one of the cell sizes is smaller than the horizon (T^-models), the 
temperatures 5T/T are almost independent of this coordinate. For instance if Rz <^1, the values 
of 6T/T are almost independent of the z-coordinate, i.e., the values of 6T/T are symmetric about 
the plane formed by the x and y-axes. This happens because of the factor 5]^/k'^ in equation (|l]). If 
we assume a power-law power spectrum with n=l, the r.m.s. value of this factor scales as A;~^/^, so 
that most of the contribution to the sum comes from small /c-values. H Rz ^1, the term Pz/Rz in 
equation (^) will be much greater than unity when 7^ 0, so the term with Pz = will dominate 
the sum. Since this term is independent of the z-coordinate, the entire sum will be approximatelly 
independent of z. In the same way, if two cell sizes are smaller than the horizon (T^-models), the 
temperatures 6T/T are aproximately independent of these coordinates. For instance if Rx,Ry <^1, 
the values of 6T/T are almost independent of both x and y, i.e., the values of 5T/T are almost 
constant along rings around the z-axis. 

The results above remain valid for a much broader range of n- values (actually, n < 3). Thus, 
the following analysis is applicable in any other large scale CMB experiment as well as one-degree 
CMB experiments. Although the existence of these symmetry patterns in the large scale fluctuations 
6T/T do not depend on the assumptions of gaussian statistics and absence of correlation between 
multipoles (see S93), we use both of these standard assumptions in this analysis. 

The analysis of T^- and T^-models is not an easy task, since there are infinitely many combi- 
nations of different cell sizes and cell orientations. In order to keep our analysis simple, we wish 
to adopt a statistic that is independent of cell orientation. In addition, we want a statistic that is 
precisely sensitive to the type of symmetries that small universes produce, so that it can rule out 
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as many false models as possible. Finally, we would like to have a statistic that is easy to compute 
and that produces results that are easy to interpret. Having these criteria in mind, we choose a 
statistic in which we calculate the function S{ni) defined by 

where Npi^ is the number of pixels that remain in the map after the Galaxy cut have taken place, 
iiij denotes the reflection of in the plane whose normal is fij, i.e., 

nij = nj - 2(nj • nj)ni (8) 

and cr(hj) and o"(njj) arc the r.m.s. errors associated with the pixels in the directions hj and hij. 
S(hi) is a measure of how much reflection symmetry there is in the mirror plane perpendicular to 
rij. The more perfect the symmetry is, the smaller -S'(nj) will be. When we calculate S{ni) for all 
6144 pixels at the positions ftj, we obtain a sky map that we refer to as an S'-map. This sky map 
is a useful visualization tool and gives intuitive understanding of how the statistic S{hi) works. 

In order to better understand S{ni), we first consider the simple model of a T^-universe 
with Rz ^1. For this specific model, the values of 5T/T are almost independent of the z- 
coordinate, so we have almost perfect mirror symmetry about the ccy-plane or, in spherical co- 
ordinates, 6T/T{9,(p) ~ 5T/T{'K — 0,(j)). When rij points in the direction of the smallest cell size 
(i.e., z-direction) , we have S{ni) ^ 1; otherwise, -S'(ni) > 1. An S-map for a T^-model {Rx,Ry,Rz) 
= (3,3,0.3) can be seen in Figure IC (lower left). Notice in this plot that the direction in which 
the cell is smallest can be easily identified by two "dark spots" at lij z and ftj — z. For 
r^-models, the only difference will be that in the place of the two "dark spots" , we have a "dark 
ring" structure in the plane formed by the two small directions. See Figure ID (lower right), an 
S-raap of the T^-model (E^, Ry, R^) = (0.3,0.3,3). 

From these two iS-maps, we can infer two important properties: first, the direction in which 
the S-map takes its minimum value, denoted 6*0, is the direction in which the universe is small. For 
a large universe such as (3,3,3), the So-directions obtained from different realizations are randomly 
distributed in the sky. Secondly, the distribution of So-values changes with the cell size, i.e., as the 
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universe becomes smaller, the values of decrease. From the definition of the S-map, it is easy 
to see that the value of So is independent of the cell orientation. In other words, if we rotate the 
cell, we will just be rotating the S-map, leaving its minimum value So unchanged. 

A value of 5*0 from a particular realization of a stochastic cosmological perturbation differs 
from the expectation value (So) due to cosmic variance. This comes from the non-symmetric part 
of 6T/T fluctuations produced by perturbation modes with Px + Py + Pz 7^ 0. For these modes, 
the main contribution to 5'(nj) in ^ is made by the terms in the sum for which and hij are 
widely separated, so that we can neglect their cross-correlation. Since S ^ \ (t") / ' 'where 

\ / ns 

{[ — )) is the r.m.s. value of the non-symmetric part of 6T/T and a is the r.m.s. noise, 

\ I ns 

we have that the cosmic variance is AS" = J (S^) — (S)'^ ^ -ta/tttt ( (^) ) ~ 0.25*. Here 



ns 



/ ~ 15 is the characteristic multipole for COBE data (the inverse angular correlation radius) and 

//.^n2\ 1/2 

\ ( ^ ) / C70 . We shall confirm this rough estimate in more details below; see the behavior 

\ I ns 

of curves for the cumulative probability distribution of So in Figure 2. 

In summary, our statistic 5*0 has all the properties that we desire: it quantifies the "smallness" 
of a sky map in a single number, it is independent of the cell orientation, and it is easy to work 
with and to interpret. 

From here on, we will present our results in terms of the cell sizes Rx^ Ry and Rz, usually sorted 
as Rx ^ Ry < Rz- We remind the reader that the results are identical for all six permutations of 
Rx , Ry and R^ . 



3. DATA ANALYSIS 

We rewrite the exponential in equation (|l|) as 

g«k-r _ k • r + i sin k • r = cos(27r7) + i sin(27r7) , (9) 

where k is given by (^, r is the vector with lenght Rh = 2cHq ^ and 7 = (^j^x + + jf-z^ . If 
the density fluctuation in Fourier space (^k has random phases, we have 5k = -^(51 + ^52)1 so that 
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(I'^kP) = (|5i+*52p-^^) = 2A^2, where is a constant and gi and g2 are two independent Gaussian 
random variables with zero mean and unit variance. Assuming a power law power spectrum with 
shape P{k) = (|(5kP) = Ak^, where A is the amplitude of scalar perturbations and n the spectral 
index, we have N = y^^"''^, so that the r.m.s. of the term 5\^/k'^ in (|l|) is given by 



n — 4 



^2 oc a 4 , (10) 

where a = (-^^ + (^) + {jt) ^ Substituting and (10) into (||), we can construct 
simulated skies by calculating 

ST 4 

— {e,<)))(x [c/icos(27r7) +fi(2sin(27r7)]a^. (11) 

Px,Py,Pz 

Since the cubic case has already been ruled out as an interesting cosmological model (see 
e.g. dOCS95), we restrict our analysis here to the and cases for n=l. This is a two parameter 
family of models specified by and Ry, with R^ = oo. For numerical convenience, we set Rz='i 
instead, as this is found to give virtually the same results as i?^ = cxd. We adopt n=l, as we found 
that "small universe" models with different n-values are even more inconsistent with the observed 
data. 

The large scale fluctuations observed on the celestial sphere by a CMB experiment can be 



modeled as being the fluctuations given in (11) multiplied by an experimental beam function 



^-{RHBk^?/2^ (12) 

where k±_ is the length of the k-component perpendicular to the line of sight and 9 is the width of 
the Gaussian beam given by 6* = FWHM/\/81n2 « 0.43 FWHM, where FWHM is the full width 
of the beam at its half maximum. We make the approximation that the sky area covered by the 
beam is flat (this is equivalent to smoothing in the plane perpendicular to the line of sight r). Since 
A;_i_ = |r X k|, we have that /c^ = fc^ — (r • k)^. We use FWHM = 7° in our simulations, which is the 
FWHM of the COBE/V)MR beam. 

In the real sky map, we do not have complete sky coverage. Because of the uncertainty in 
Galaxy emission, we are forced to remove all pixels less than 20° below and above the Galaxy plane. 
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which represents a loss of almost 34% of all pixels. However, performing Monte Carlo simulations 
with and without the Galaxy cut, we find that the Galaxy cut does not change the final distribution 
of 5*0 much. Due to the smaller data sample, the Galaxy cut weakens the lower limit on the cell 
size slightly (see e.g. Scott et al. 1994). 

We model the noise rii at each pixel i as independent Gaussian random variables with mean 
(rij) =0 and variance {niUj) = (Jij^ij (Lineweaver et al. 1994), and add it to the temperature values 
given by ([Tl|). The level of noise in the DMR maps is a source of serious concern in our analysis: 
high levels of noise can make it impossible to discriminate between the different topological models. 
In order to reduce the noise and increase the signal-to-noise ratio in the simulated skies and real 
data, we smooth both once more by 7° before calculating S-map, which corresponds to a total 
smearing of V(7°)^ + (7°)^ « 10°. 

We generate our simulated skies as standard DMR maps with 6144 pixels for n=l, with a 
Galaxy cut of 20°, FWHM = 7°, the monopole and dipole removed, add noise and normalize to 
(T70 = 34.98^K (the r.m.s. value at 7° extracted from our DMR map, a 4 year combined 53 plus 
90 GHz map with monopole and dipole removed). Fixing a cell size, we construct a simulated sky 
according to ([TTI), we smooth this once more by 7° and use the statistic defined in (0) to obtain 
an S'-map from which we extract its minimum value So. Repeating this procedure 1000 times, 
we obtain the probability distribution of 5*0 for that fixed cell size. When we repeat this same 
procedure for different cell sizes, we are able to construct Figure 2. 

In Figure 2A (upper plot), we show the cumulative probability distribution of So obtained from 
the Monte Carlo simulations for the cell sizes {Rx, Ry, Rz) = (0.5,0.5,3), (0.6,0.6,3), (0.7,0.7,3) and 
(3,3,3). The horizontal lines indicate the confidence levels of 95%, 90% and 68% (from top to 
bottom). Comparing these curves with the value 5,^^^^ = 2.59 (represented in the plot by the 
vertical straight line), where S^^'^^ is the 5o value extracted from our data set, we conclude that 
T^-models with smallest cell sizes Rx,Ry ^0.5 can be ruled out at 95% confidence. As the second 
cell size Ry is increased, the curves shift to the left of the T^-models and we can rule out T^-models 
for Rx ^0.5 at a similar confidence level, see Figure 2B (lower plot). In this plot, we show the 
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cumulative probability distribution of So obtained from Monte Carlo simulations for the cell sizes 
{R^,Ry,R,) = (0.5,3,3), (0.6,3,3), (0.7,3,3) and (3,3,3). 

A more complete picture of the cell size limits is obtained when we construct a two-dimensional 
grid for different values of the cell sizes {Rx,Ry,Rz) with R^ = 3.0 and 0.2 < Rx,Ry < 3.0 (see 
Figure 3). The thin-shaded, thick-shaded and grey regions correspond, respectively, to the models 
ruled out at 68%, 90% and 95% confidence. Notice in this plot that all contours are almost L- 
shaped, which means that, to a good approximation, the level in which a model {Rx,Ry) is ruled 
out depends only on the smallest cell size, Rmin = {Rx,Ry}- We see that Rmin ^0.5 at 95% 
confidence. 

4. CONCLUSIONS 

We have shown that the COBE /DMR maps have the ability to discriminate and rule out 
and topological models. We have presented a new statistic to study these anisotropic models 
which quantifies the "smallness" of a sky map in a single number. So, which is independent of the 
cell orientation, is easy to work with and is easy to interpret. 

Prom the COBE/DMR data, we obtain a lower limit for T^- and T^-models of R^ >0.5, 
which corresponds to a cell size with smallest dimension of L=3000/i~^Mpc. This limit is at 95% 
confidence and assumes n=l. Since the topology is interesting only if the cell size is considerably 
smaller than the horizon, so that it can (at least in principle) be directly observed, these models 
lose most of their appeal. Since the cubic case has already been ruled out as an interesting 
cosmological model (see e.g. dOCS95), and T^- and T^-models for small cell sizes are ruled out, 
this means that all toroidal models (cubes and rectangles) are ruled out as interesting cosmological 
models. 
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FIGURE CAPTIONS 

Figure 1: Simulated sky maps for the T^- and T^-models and their /S-maps. (A) r^-model with 

dimensions {Rx, Ry, Rz) = (3,3,0.3); (B) T^-model with dimensions (R^, Ry, Rz) = (0.3,0.3,3); (C) 
S-map of the T^-model shown in (A); (D) S-map of the T^-model shown in (B). 

Figure 2: Cumulative probability distribution of So for T^- and T^-models obtained from Monte 
Carlo simulations. (A, upper plot) Simulations for T^-universes with dimensions {Rx,Ry,Rz) = 
(0.5,0.5,3) or dot-dashed line, (0.6,0.6,3) or dashed line, and (0.7,0.7,3) or dotted line. (B, lower 
plot) Simulations for T^-universes with dimensions {Rx, Ry, Rz) = (0.5,3,3) or dot-dashed line, 
(0.6,3,3) or dashed line, and (0.7,3,3) or dotted line. In both pictures the model {Rx, Ry, Rz) = 
(3,3,3) is represented by a solid line, S^^^ = 2.59 (vertical straight line) and the horizontal solid 
lines indicate the confidence levels of 95%, 90% and 68% (from top to bottom). 

Figure 3: Grid of cumulative probability distributions of 5*0 for T^- and T^-models obtained from 
Monte Carlo simulations. The thin-shaded, thick-shaded and grey regions correspond, respectively, 
to the models ruled out at 68%, 90% and 95% confidence. 
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